The Damping of the Bose-Condensate Oscillations in a Trap at Zero Temperature 
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We discuss an existence of the damping for the radial condensate oscillations in a cylindric trap 
at zero temperature. The damping is a result of the parametric resonance leading to energy transfer 
from the coherent condensate oscillations to the longitudinal sound waves within a finite frequency 
interval. The parametric resonance is due to the oscillations of the sound velocity. The triggering 
amplitudes at zero temperature are associated with the zero-point oscillations 



INTRODUCTION 

The damping of the oscillations of a trapped Bose con- 
densate isolated from environment is one of the most in- 
teresting problems in the physics of the Bose-Einstein 
condensation. So far, experimental [1-4] and theoret- 
ical [5-7] investigations have been reduced to studying 
the damping due to the interaction of oscillations with 
thermal normal excitations. The ensemble of such ex- 
citations in essence plays a role of the heat bath. In 
all cases the relatively high temperatures T ^ fiwo have 
been considered, being the frequency of a parabolic 
trap. However, the principal question about the origin of 
the irreversible damping in the oscillating isolated system 
at r = remains open. 

In the present work it is shown that such damping at 
zero temperature, at least at some definite conditions, 
really exists. We consider the radial oscillations of the 
condensate in an elongated trap with the cylindric sym- 
metry, which are induced by the rapid decrease of the 
frequency Lo{t) of the transverse isotropic potential from 
value ojQ to lji. Such statement of the problem has a 
number of advantages. As is found in [8], there exists 
an exact scaling solution of the nonlinear Schrodinger 
equation for a field operator and its quasiclassical ana- 
log, Gross-Pitaevskii equation, for an arbitrary behavior 
Loit) in the case of a gas in the isotropic two-dimensional 
parabolic potential. It is essential that this solution holds 
for the quasi-2D case of an elongated trap with the cylin- 
dric symmetry. In particular, a ratio between the interac- 
tion energy of particles and their kinetic energy remains 
constant. Thus, if the Thomas-Fermi approximation is 
applicable in the initial static trap, it remains valid at all 
stages of the gas evolution. 

The scaling solution describes the space-time evolution 
of the condensate, relying only on the solution in the ini- 
tial static potential of frequency loq. After the parameter 
of the potential takes the value wi, the condensate sets 
into the stable state accomplishing with the radial oscil- 
lations of frequency 2cji. These oscillations are accom- 
panied by the vibrations of the condensate density and, 
therefore, by the vibrations of the sound velocity c. 

As will be shown below, for a gas in which the sound 
velocity oscillates, the phenomenon of parametric reso- 
nance appears. The essence of the phenomenon is that 



the resonance between the radial condensate vibrations 
and the longitudinal wave propagating in the direction of 
the axial z axis produces an exponential increment for the 
amplitude of the wave. At zero temperature the trigger- 
ing amplitudes are formed, in essence, by zero-point os- 
cillations. In the Thomas- Fermi regime, when the chem- 
ical potential /i 3> hw^, the sound waves of wavelength 
A ~ c/wi amplifies. As a result, the dynamic energy of 
the coherently oscillating condensate reduces, converting 
into the energy of the longitudinal sound vibrations. For 
the analysis of the initial stage of generation at T = 0, it 
is of interest to find the quantum-mechanical description 
of the process. Accordingly, in the last section the rate 
is found for the elementary process of creating a pair of 
longitudinal phonons in a periodic field of the radially 
oscillating Bose-condensate. 



THE BOSE-CONDENSATE OSCILLATIONS IN A 
TRAP 



Let us consider a Bose-gas in the elongated cylindric 
trap with the symmetrical parabolic potential and lon- 
gitudinal size L ^ Rq. Neglecting edge effects, we can 
represent the general equation for the Heisenberg field 
operator 5* (r, z, t) in the cyhndric coordinate system as 



ih 



5^- 

dt 



2m 



2m dz^ 



f7o*+**. (1) 



Here Uo — Anah^/m where a is the scattering length. 
The single simplification in (|l]) is an assumption on the 
local character of the interparticle interaction. 

A unique property of the above equation, contain- 
ing parabolic 2D potential of frequency uj{t) with an 
arbitrary dependence on the time, is an existence of 
the scaling transformation. This transformation reduces 
equation (|^) to the form in which the frequency of the 
parabolic potential is constant and equal to the initial 
one LjQ. 

Following [8], let us introduce new scale p = r/b{t) in 
the transverse direction and simultaneously define new 
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time variable r(t). Representing the field operator as 

^ (r, z, t) = [p, z, t) exp [^$ (r, i)] (2) 

we insert this expression into (|^) . In terms of new vari- 
ables the equation for operator x can be reduced to (see 
Appendix) 

dx . 2 1 2 2 1 - 

+ C/oX-XX-5^w|^g, (3) 

if the phase is expressed via the parameters of trans- 
formation h (t) and T{t) by the relation 



mr^ db IJ,T{t) 



2hb dt h 
and these parameters satisfy the equations 



II = const, (4) 



dt 



dt 



(5) 



(6) 



Here luq = uj{—oo) and Eq. (^) has the following initial 
conditions b{t ^ — oo) = 1. In the case of the prob- 
lem uniform in the z direction the equation in variables 
p, T reduces practically to the equations for the static 
2D parabolic potential of frequency ujq. Assuming rapid 
switching at the initial time moment for the trap fre- 
quency from ujq to coi, we find for the solution of 
Eq. ® 

W = ^ + 1) (1 -9Cos2LUit), t>0 (7) 

where (3 = ujq/uji > 1 and g = (/J^ - l) / (/J^ + l). Thus 
the parameter b{t) oscillates at frequency 2uji within the 
interval from 1 to /3. 

Substituting this solution into (|^) and integrating, we 
obtain 



term on the r.h.s. of describing free motion of the gas 
along the cylindric axis. 

Considering the ground state at T = 0, we can replace 
as usually the operator x in (|3|) with the macroscopic 
wavefunction of condensate xo- Taking into account that 
the symmetry of the problem determines an indepen- 
dence of Xo on z, we arrive at the Gross-Pitaevskii equa- 
tion, e.g., [10] 



.,dxo ,^^ 2,1 2 2 



2m 



fAxo + UoXoXoXo- (9) 



From the general representation with one can 
conclude that /i is a chemical potential corresponding to 
the static trap of frequency uq. In addition, the left-hand 
side of Eq. vanishes and xo in variables p and r is a 
stationary real solution of the equation 

i-^K + l^^op' -^^]xo + Uoxl = o (10) 

Let us restrict ourselves by the case when the inequality 
holds for 



/i > Hojo 



(11) 



and therefore the Thomas- Fermi approximation realizes. 
Then one can neglect the "kinetic energy", namely, the 
first term in Eq. (|l0|), and obtain the known solution 



Xo 



1/2 



^1 (1-7^ 



P 



1/2 



where 



Rq = J 2/i/ mujl 



(12) 



(13) 



The number of particles and the energy corresponding to 
this wavefunction are equal to 



iVo 



4a \ fujjr) 



(14) 



The energy of the ground state in the trap of frequency 
uji equals 



loqt {t) = ujit + arctan[/3tan (wit)] 

— arctan[tan (cjii)], t>0 (8) 

At the stages of small expansion the time r varies at the 
rate closed to the laboratory one t. For the stages of 
the maximum expansion, the time r in the frame moving 
together with the gas varies slower. 

In the general case the nonlinear operator Schrodinger 
equation (||) describes the non-ideal Bose-gas in the con- 
comitant system of coordinates (p, r) in which the gas 
is in the static potential. A single explicit manifestation 
of the radial condensate oscillations contains in the last 



El = -MiA^o 



where 



Pi = P- 



I — . 

UjQ 



(15) 



(16) 



It is worthwhile that solution ( |12| ) is valid at the large 
negative times when uj(t) = loq as well as at the large 
positive times when u;(t) — lui. This means that, on one 
hand, the function xo is the wavefunction of the ground 
state of a gas in the trap of frequency loq and, on the other 
hand, the same function describes the condensate state at 
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the large positive times. The state Xo is static only in the 
variables p and r but in the laboratory frame it oscillates 
with the frequency 2uji as follows from (Q). To describe 
the space-time evolution of the condensate wavefunction 
in the variables r and t, it is sufficient to substitute ( p^ ) 
into the quasiclassical analog of expression (||) and use 
(D, (0), and (D. 

EXCITATIONS OF THE BOSE-CONDENSATE 

With respect to the trap with the parabolic potential 
of the finite frequency oji , the oscillating state xo is an 
excited state. Naturally, the question arises whether this 
state is decaying and if it does, what mechanism is. Let 
us consider the excited states of the system at the back- 
ground of the coherently oscillating condensate. These 
oscillations can be found as oscillations of the classical 
field of the condensate, see, e.g., [11]. For this purpose, 
let us introduce function X = Xo + x' instead of operator 
X in Eq. (|^) and linearize it in x' 



Inserting ( pO| ) into Eq. ( |l9| ) for x' and x'*j we find 



ih 



dF b^h^P 



dr 



OT 



2m 
2G- 



f, 

2m 



F. 



Let us go over in these equations from variable t to the 
laboratory time t, using relation (^. Then the equations 
can be represented as a second-order equation 



d^F 



2m 



2G_ 

WW) 



2m 



F = 0, 



(21) 



Let us average Eq. (^ ) over p and introduce notation 
(F) = F. The average of a product GF, involving that 
G is a positively determined function, can be represented 
as 



(GF) = GF 



As a result, we obtain 



_dx' 



1 



2 

G(x' 



mujlp^ + G - p\x' 



2m 9z2 • 



Here 



G = Uoxl 



(17) 



(18) 



We are interested in the long wavelength longitudinal ex- 
citations. For fcz = fc — > 0, as is clear from the physical 
reasons, the lowest branch of the longitudinal excitations 
is associated with the uniform shift along the z axis with 
the radial distribution of the density determined by the 
function xo- Hence this branch is gapless. 

For the small finite k, the transverse distribution of 
the density in such wave changes weakly. This is clearly 
seen, e.g, from the results obtained in [12,13] for a static 
elongated trap. Taking this into account and employ- 
ing Eq. ( [l0| ) and definition (|l8|), it is easy to find that 
the first term on the r.h.s. of ( p7| ) is obviously small 
compared with the term Gx' ■ Then Eq. (17) simplifies 
significantly 



^n^ = G{x' + x'n 



fi2 «2 / 



2m 



dz 



(19) 



We seek for the solution of the equation for the longitu- 
dinal excitation in the form, see, e.g., [11] 



X'(p,z,r) = u(p,r)e'*=^-w*(p,T)e 



~ikz 



(20) 



Let us introduce the notations 
1 



/ 



(u + v), F 



Xo 



(u — v) 



d^F 



+ nl it)F = o, 



2m 



2G 



62 {t) 2m 



(22) 
(23) 



For the case of the static potential b — 1, flk is nothing 
else than the frequency of the Bogoliubov spectrum for 
the longitudinal excitations, determined for some average 
condensate density. The sound velocity equals 



G 

Co = \/ — 
V m 



(24) 



In the problem which we consider, the scaling parame- 
ter b oscillates in time. As follows from (p3|), this results 
in the oscillation of the sound velocity. As will be seen 
below, we are interested in the longitudinal phonons of 
frequencies Qk — (^i- Owing to inequality ( pT] ) this is the 
purely acoustic region. Accordingly, 



(25) 



Using explicit form of the dependence b{t) (0), we 
have for the coefficient ^1{t) in Eq. (I2F 



it) 



{I — g cos 2LJit) 



where 



2 (cfc)' 
13^ + 1 



(26) 



(27) 
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Thus, Eq. ( p2D is reduced to the known Hill equation. 
For the relatively weak variation of the trap frequency 
when (7~/3— !<<!, we arrive at the Mathieu equation 



+ ujI{1 + g cos 2ujit) F = 0. 



(28) 



This equation, as well as a more general Hill equation, de- 
termines the parametric resonance connecting the coher- 
ent transverse condensate oscillations at frequency 2wi 
with longitudinal phonons with the frequencies close to 
uji. From the physical point of view the parametric res- 
onance arises due to periodic variation of the condensate 
density and, therefore, of the sound velocity (G£^ 



THE PARAMETRIC RESONANCE 

To clarify the conditions for appearing the paramet- 
ric resonance, we restrict ourselves by the case g <C 1 
and employ the standard algorithm for solving Mathieu 
equation, see e. g. [9]. We seek for a solution as 



F ^ A{t) cosLUit + B (t) sin Wit. 



(29) 



The coefficients A{t) and B{t) are assumed to be the 
slowly varying functions of the time t and on the substi- 
tution of (|9|) into the Mathieu equation (^ ) one should 
retain only the terms zero and linear in g, neglecting sec- 
ond derivatives of A{t) and B{t). As a result, we arrive at 
a set of two algebraic linear equations with the constant 
coefficients. One of the solutions grows exponentially. 



A{t),B{t) 




(30) 



provided that the frequency lies within the narrow band 
near uji with the width 



6uj = t^^ = ^5^1- 



(31) 



Within this band a growth of the amplitude of the sound 
waves takes place and therefore the energy converts into 
the longitudinal modes. As the result the energy of the 
coherent transverse condensate oscillations decreases. It 
is necessary to have at least one mode of the longitudinal 
sound excitations within the energy band of the para- 
metric resonance (^). The spacing between the neigh- 
bor modes equals Alu — 2nc/L. Since Rq — 2c/ wi, this 
condition holds for if there is an inequality 



L > 2TTRo/g. 



(32) 



Within the narrow band around lui a square of the 
wavefunction amplitudes ( pO| ) increases as exp 2^kt- At 
T = the initial energy equals the energy of a mode of 



zero-point oscillations. Thus the total energy transferred 
into the longitudinal modes due to parametric resonance 
equals 



E{t) 



^ 2 

fe 



(33) 



Let us suppose that a sufficient number of modes lies 
within the band (31) and a sum in ( p3| ) can be replaced 
with an integral 



Jit) 



dx 



gLhhJi 



IGttc 
t 



J{t), 



(34) 



exp 



Let us compare the energy (|^ with the energy of the 
coherent condensate oscillations. For an instantaneous 
transition of the trap frequency from ujq to toi , the wave- 
function xo does not change its configuration. Hence the 
energy of the condensate falls from the initial magnitude 
Eq (see (|l^) to the magnitude 



1 En 



E'A{% 



The energy of the vibrational motion E^ which the con- 
densate possesses is equal to a difference between this 
energy and the energy of the ground condensate state 
( [Tsl ) corresponding to frequency 



E,^E' ^Ei = -fiNog^ 



(35) 



As a result, the damping of the condensate oscillations is 
characterized by a ratio 



E{t) SaujQ 
Airgc 



Er 



Jit) 



a 
gRo 



J{t). (36) 



For the small times J(t) 
duces to 



t/ti, and the ratio (36) re- 



3 aujQ 



hujQ 



t. 



(37) 



Ec 16 c 

For the times t 3> ii, the damping grows exponentially 



m 

Er 



2tt gRo 



2TTtl 



1/2 



exp 



(38) 



The factor in front of the exponential is much smaller 
compared with unity for the realistic values of 5. Hence, 
within the logarithmic accuracy the typical time for de- 
caying the condensate oscillations is equal to 




(39) 
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This time is large compared with the period of conden- 
sate oscillations not only due to the smallness of g but 
also due to the large factor in the argument of loga- 
rithm. Thus the damping of the condensate oscillations 
due to the mechanism considered is a rather slow process 
though appearing inevitably at T = under conditions 
concerned. 

The parametric resonance leads to reducing the num- 
ber of particles in the condensate. In principle, this could 
require to solve a self-consistent problem. However, from 
the general considerations it is clear that the escape of 
particles from the condensate can be neglected if the vi- 
brational energy per particle g^/i (see ( |35| ) ) is small com- 
pared with the temperature of Bose-condensation Tc- 



THE QUANTUM CALCULATION OF THE 
PHONON GENERATION AT THE INITIAL 
STAGE 



The classical parametric resonance develops from the 
moment t = only in the case when nonzero amplitude 
of the resonant modes within the interval 6uj (^Tj) exists 
from the very beginning. At the finite temperature these 
are amplitudes of thermal oscillations. At T = the ini- 
tial amplitudes should be connected with the zero-point 
oscillations of the longitudinal modes. This is explicitly 
taken into account in the preceding section. It is of in- 
terest to trace the initial stage of the phonon generation, 
remaining within the framework of the purely quantum 
mechanical consideration. For this purpose, we should 
return to the operator description of excitations in the 
system. 

The general equation for operator x' lias the form anal- 
ogous to that of Eq. (UTv) 



^f^^-h,x' + G{x' + x'^)-^b^r)% (40) 



r 2m oz'^ 



(41) 



We go over from variable r to t (see (^)) and single 
out the term associated with the coherent condensate 
oscillations in the explicit form 



dx' 



zft^ = hx' + [h,x' + G{x' + {b-^ [r) - l) 



(42) 



2m dz^ 



h = hn 



Assuming that (6^^ (r) — l) << 1, we treat the last 
term on the r.h.s. as a perturbation. From the stan- 
dard commutation rule for the field operators \l/(r, z,i) 
and representation (0) it directly follows the commuta- 



tion rule for operators x' 

[X'{p,z,t),x'^{p\z',t)] 

= b^6 (r - r') S{z-z')^6{p- p') S {z - z') 

Using this relation, one can readily write the perturba- 
tion hamiltonian governing the behavior of the last term 
on the r.h.s. of Eq. (42) 



H'={b-^ (r)-l) 



d pdz 



1 



X'+ (/io + G)x' + ^G(x'x' + x'+x'+) 



(43) 



Consider at first the excitations in the static case when 
b — 1. In the secondary quantization the operator x' can 
be written using the eigenfunctions of hamiltonian h as 



X' = ^ankXnk, Xnk = -j=LPn{p)-, (44) 

nk 

holfinip) = Enipnip) 

Inserting (^J) into Eq. (^), multiplying both sides of 
the equation by x*^., and integrating over p and z, we 
find 



2m 



where 



G - ' "'"2 



d'p^i{p)G{p) 



Equation ( [is] ) is approximate since we have omitted the 
terms nondiagonal in n and proportional to Gnn' ■ The in- 
volvement of these terms leads to the quantitative correc- 
tions alone. Employing the standard Bogoliubov trans- 
formation 



a„fc — Unkbnk ~ Vnkb^_,., 

and obvious condition 



(46) 



ih- 



dbn 



dt 



£nkbr, 



for independent collective excitations, we find the spec- 
trum of excitations and the values of coefficients Unk, Vnk- 
For the reasons mentioned in the preceding section, we 
consider only the phonon generation corresponding only 
to the lowest branch of excitations n = 0. For Enk in this 
case, we find the expression coinciding with ( |23| ) at & = 1 
if Goo is understood as G. Now let us find the phonon 
generation in first approximation in H' (^). For this 
purpose, using representation (|4^), we perform transfor- 
mation (|4q), considering only the terms with n=:0. Since 
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we consider the case T=0, the corresponding generation 
proves to be associated only with the terms having a 
product bQi^b^_j. in the hamiltonian, i.e., with the produc- 
tion of two phonons with the opposite momenta. Corre- 
spondingly wc find directly for the transition amplitudes 



Ak-k = -^gG cos (2wii) {uok - vok) 

X (fc, ~k 



bokbtk 



0) (47) 



Thus, at the initial stage of relaxation the probability of 
generating phonon pairs per unit time is equal to 

Wk,-k = ^ [gGiuok - vokf] ' S{2hui - 2eQk). (48) 

Using the known relations for the coefficients of the Bo- 
goliubov transformation, e.g., [10], we have for the region 
of sound excitations (see (|27|)) 

EOfe 

Hence we readily find the power of the energy lost by the 
oscillating condensate 

L f , Lh 



(uok - VQk ) 



eok — fiiok 



E : 



277 fic 



de2eWk 



64c 



(49) 



This result coincides with that ( |37| ) (taking into account 
( ^ ) and (p^) obtained in the preceding section for the 
parametric resonance at the initial times f <C ti. The 
quantum mechanical calculation confirms our notion on 
generating the longitudinal waves at zero temperature as 
an enhancement of zero-point oscillations due to para- 
metric resonance. 



THE PARAMETRIC RESONANCE AT THE 
MULTIPLE FREQUENCIES 



As is known, the Mathieu equation ( p8[ ) results in the 
parametric resonance at the frequencies close to the mul- 
tiple frequencies poJi. However, considering the phonon 
generation at these frequencies, we must return to the 
general equation ( |2^ ) even in the case 5 <C 1. The anal- 
ysis of the equation shows that the parametric resonance 
is absent for the even multiple frequencies. The phonon 
generation can occur only near the multiple frequencies 
puji with odd p. As p grows, the growth decrement E(t) 
and the width of the resonant interval decrease strongly 
at any value of g. This can easily be traced, analyzing the 
nearest multiple resonance of p = 3. Such unusual result 
is due to specific behavior of the oscillating coefficient in 
Eq. (G3) varying periodically in time. 



CONCLUSIVE REMARKS 

Thus, analyzing the coherent radial oscillations of the 
condensate in a cylindric trap, we provide the evidence 



for the existence of damping of these oscillations at zero 
temperature. The damping is a consequence of the para- 
metric resonance caused by the sound velocity oscilla- 
tions. As a result the energy transfer from radial con- 
densate oscillations to longitudinal phonon modes is re- 
alized. 

At the finite temperature hwi < T << /i the same 
resonance amplification takes place but now the initial 
number of phonons with the energy ftiOi equals T/hwi 
and the initial energy of the resonance modes equals hwi ■ 
{T/hui). For this case, in Eq. ( |3^ ) the quantity huJi 
should be replaced with T. Such replacement increases 
insignificantly the argument in logarithm (|39|). 

Let us estimate the characteristic damping time for 
quantum gases Rb and Na. We consider an elongated 
cylindric trap with j/q = 400Hz (Ru) or i'q — lO'^Hz 
{Na), L = 2 ■ 10^ and assume that the dimensionless pa- 
rameter g w 0.15. Then from Eq. ( |l4|) we find /i ~ 90nK 
for the Rb gas. Accordingly, the Thomas-Fermi radius 
equals R « 1.7 • 10~'*cr7i. For the parameter ti (^l]), tak- 
ing into account the relation uji — (1 — g) ■ coq, we have 
ti Ri 6.2 • lO^'^s. The direct calculation estimates the log- 
arithmic factor in Eq. ( ^ ) as ~ 10. Thus, evaluating the 
damping time, we find < OAs. Since it is smaller than 
the ordinary lifetime of the system, this time is realistic 
for revealing the damping. 

The similar results can be obtained for sodium. In this 
case the estimates give fj, w lAOnK, R w 1.6 • 10"'*cm, 
and ti « 2.5 • 10~'^s. As a result, t* ~ 0.03s. It is impor- 
tant that in both cases the condition (^2|) proves to be 
fulfilled. At the same time gfi <^ Tc, and this means that 
the depletion of the condensate is quite small after the 
complete damping. 

This work is supported by the Russian Foundation for 
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tion for Scientific Research (NWO). 



[1] 
[2] 

[3] 

[4] 

[5] 
[6] 

[7] 

[8] 

[9] 



D. S.Jin, J.R Ensher, M.R. Matthews, C.E. Wieman and 

E. A. Cornell, Phys. Rev. Lett. 77,420 (1996) 

M.O. Mewes, M.R. Andrews, N.J. van Druten, D.M. 
Kurn, D.S. Durfee, C.G. Townsend and W. Ketterle 
Phys. Rev. Lett. 77, 988 (1996) 

D. S.Jin, M.R. Matthews, J.R Ensher, C.E. Wieman and 

E. A. Cornell, Phys. Rev. Lett. 78, 764 (1996) 

D.M. Stamper-Kurn, H.-J. Miesner, S. Inouye, M.R. An- 
drews and W. Ketterle, Phys. Rev. Lett. 81, 500 (1998) 
W.V. Liu, Phys. Rev. Lett. 79, 4056 (1997) 
L.P. Pitaevskii and S. Stringari, Phys. Lett. A235, 398 
(1997) 

P.O. Fedichev, G.V. Shlyapnikov and J.T.M. Walraven, 
Phys. Rev. Lett. 80, 2269 (1998) 

Yu. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 
Rev. A 54, R1753 (1996) 

L.D. Landau and E.M. Lifshitz, Mechanics, (Pergamon 
Press, Oxford, 1980) 



7 



[10] E.M. Lifshitz and L.P. Pitacvskii, Statistical mechanics, [12] E.Zarcmba, Phys. Rev. A 57, 518 (1998) 

Part 2, (Pcrgamon Press, Oxford, 1988) [13] S. Stringari, Phys. Rev. A 58, 2385 (1998) 

[11] Fr. Dalforo, St. Giorgiiii, L.P. Pitaevskii and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999) 



